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ABSTRACT 

The accurate determination of cluster total mass is crucial for their use as probes of cosmology. 
Recently, the Sunyaev-Zel'dovich effect (SZE) has been exploited in surveys to find galaxy clusters, 
but X-ray or lensing follow-up observations, or empirically determined scaling relations between SZE 
flirx and total mass, have been required to estimate their masses. Here, we demonstrate a new 
method of mass determination from SZE observations, applicable in the absence of X-ray or lensing 
data. This method relies on the virial relation and a minimal set of assumptions, following an approach 
analogous to that used for stellar structure. By exploiting the virial relation, we implicitly incorporate 
an additional constraint from thermodynamics that is not used in deriving the equation of hydrostatic 
equilibrium. This allows us to relate cluster total mass directly to the robustly-determined quantity, 
the integrated SZE flux. 

Subject headings: cosmology: observations — dark matter — galaxies: clusters: general — galaxies: 
clusters: individual (A1835, A1914, CL J1226.9+3332) 



1. INTRODUCTION 

Clusters of galaxies are thought to be the largest 
gravitationally-bound objects in the universe and there- 
fore good tracers of cosmology. Ongoing cluster 
surveys, such as th ose with the A tacama Cosmol- 
ogy Telescope (ACT: lKosowsk^[200l 
^ninll. the South P ole Telescope (SPT^ 
Plagge et a l.' '2011 



Menantcau ct alj 
lR"uhl ct al. 2004 

, and Planck (jRosset et al] 126101 : 

Planck Collaboration et ani2011bl lal). have the potential 
to place tight constraints on cosmological parameters 
with the clusters they discover. These surveys uti- 
lize the Sunyaev-Zel'dovich effect (SZE), which has red- 
shift independent surface brightness and arises by Comp- 
ton scattering of cosmic microwave background (CMB) 
phot ons off of the hot electron s in clusters of galax- 

y Zlel'dovich fc SunvaevI 119691 : iSunvaev fc Zel'dovichl 
However, the interpretation of cluster yields 
relies on the accurate determination of the scaling 
between integrated SZE flux and cluster total m ass 
(IMcCarthv et al.l[2003l : IMotl et al.l[2005l : WasaMM\ . 

Because the SZE intensity varies as the line-of-sight in- 
tegral of thermal electron pressure, a cluster's total, in- 
tegrated SZE flux scales linearly with the volumetric in- 
tegral of thermal pressure, which is thermal energy (see, 
e.g.. Equation H]) . To the extent that clusters are virial- 
ized and supported hydrostatically by thermal pressure, 
a cluster's total SZE flux will closely track its gravita- 
tional energy, thereby motivating SZE flux as a proxy 
for cluster mass. Approaches exploiting this expected 
tight correlation have traditionally relied on empirical 
relations between SZE flux and total mass determined 
from X-ray observations or optical lensing studies. 

We present a new approach for mass determination 
from SZE data alone, applicable in the absence of X-ray 
or lensing data, that exploits the virial relation and a 
minimal amount of simplifying assumptions about clus- 

^ Department of Physics and Astronomy, University of Penn- 
sylvania, Philadelphia, PA 19104, USA 
^ Einstein Postdoctoral Fellow 



ter astrophysics. In Section [2] we describe how SZE ob- 
servations can constrain thermal energy, and we relate 
this directly to cluster total mass via the virial theorem. 
In Section [3l we demonstrate this mass determination on 
previously published SZE data. In Section HI we assess 
how these simplifying assumptions impact this method 
and offer conclusions. 

2. THERMAL ENERGY CONSTRAINTS FROM 
OBSERVATIONS OF THE SUNYAEV-ZEL'DOVICH 
EFFECT 

The thermal SZE is a smaU (< lO^'^) distortion in 
CMB intensity caused by inverse Compton scattering 
of CMB photons by energ etic electrons in the hot in- 
tracluster medium (ICM; Zel'dovich fc SunvaevI 119691 : 
ISunvaev fc Zel'dovichl Il972t) . This spectral distortion 
can be expressed, for dimensionless frequency x = 
hv/kj^TcMB, where h is Planck's constant, v is frequency, 
fee is Boltzmann's constant, and Tqmb is the primary 
CMB temperature, as the change in intensity A/sze rel- 
ative to the primary CMB intensity normalization Iq, 
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The factor g{x,Te) in Equation [T] encapsulates the fre- 
quency dependence of the SZE intensity. For non- 
relativistic electrons, 
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(2) 



At low frequencies (< 600 GHz), relativistic corrections 
to E quation |2] are fa irly straightforward to apply (see, 
e.g., Iltoh et all |1998() . The Compton y parameter in 
Equation [1] is defined as 
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where CTt is the Thomson scattering cross-section of the 
electron, £ is the line of sight, and mgC^ is an electron's 
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rest energy, and the primary CMB intensity normaliza- 
tion is /o = 2{k^Tc„^f{hc)-^ = 2.7033 x 10^ Jy Sr-^. 
Note that we have assumed the ideal gas law (Pg = 
nefesTe) in Equation |3] to relate electron pressure Pg to 
electron number density and temperature T^. 

From Equation |3l one ean see that resolved observa- 
tions of the thermal SZE from a cluster can be used to 
constrain its electron pressure profile PcC*")- This can 
be related to the total pressure as Pgas = (1 + l/^o)^c, 
where /io = 2/(1-1- X) is the mean particle weight per 
electron and X is the mass fraction of hydrogen. While 
deep X-ray observations have shown that the distribu- 
tion of heavy elements in the ICM varies w ith radius (e.g. 
iVikhlinin et'alll2005l: [Peterson et al.ll200lD . and theoret- 
ical studies indicate that helium sedi men tation in to the 
cluster core will als o impact /ic (e.g. iMa rkevitchl 120071 : 
iPeng fc Nagail [20091 ) . we make the simplifying assump- 
tion that /Xe = 1.17. We assess the impact of this as- 
sumption in Section [H 

A common SZE observable used in SZE mass scaling 
relations is l^nt, the Compton y parameter integrated 
over some region of the sky, defined as 



y d£L. 



(4) 



Because lint is proportional to the surface brightness of 
the SZE integrated over a region of the sky, it tracks the 
integrated cluster SZE flux. For a spherically symmetric 
electron pressure profi le Pc(r), the spherically- in tegrated 
version of i^nt is (e.g.. iMroczkowski et al.1l2009[ ) 

r,ph(r) = ^ /Pe(r')47rr'2dr' (5) 

P^^{r')AW^dr' (6) 



(l + l/Me)mec2io 

2gTgth(r) 
3(1 + 1/Mc) rricC^ 



(7) 



We have used the fact in Equation [7| that the thermal 
energy within r is 



(8) 



for a monatomic, ideal gas. 

Through Equation [3 one can see that the SZE observ- 
able Kjph relates directly to the thermal energy content 
of the cluster. To the extent that a cluster is virialized 
and supported by thermal pressure, this quantity will 
closely track the gravitational energy Ug (r) via the virial 
relation, 

2Sth(r)-3P(r)y = -C/g(r). (9) 

The —iPV term, where V = Airr^ /i is the volume at 
r, accounts for the non-vanishing surface pressure and 
- as noted recently by Colin Hill - must be taken into 
account when solving for the mass. This term works 
to reduce the amount of gravitating mass required to 
hold the gas within r. Using Equation [3 which de- 
rives from statistical mechanics, to relate pressure to 
thermal energy, we note that the virial relat i on (Eq ua- 
tion [9]) is derived (see , e.g., ISchwarzschild[ ([1951 or 
[Kippenhahn &: Weigertj (|1990[ ). who derive and discuss 



the surface pressure term) from the equation of hydro- 
static equilibrium (HSE), 
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where Pga.s{i') is the gas density as a function of radius 
r, Mtotii") is the total mass within r, and G is the grav- 
itational constant. The equation of HSE is derived from 
fluid mechanics, specifically the equations of motion and 
continuity (see, e.g.. iSchwarzschild,1958; ,Sarazin.l988, ). 
Mass estimates based on HSE traditionally assume a 
spherically-symmetric mass distribution and the ideal 
gas law, and therefore adopting the virial relation is no 
more restrictive than the assumptions typically required 
for HSE mass determinations. However, Equation [5] pro- 
vides an additional, key constraint not used in pure HSE 
mass determinations, which (assuming the ideal gas law) 
only require two of these three ICM profiles: density, 
temperature, and pressure. As we discuss in iJJl the virial 
mass estimate is proportional to the square root of scalar 
changes in the pressure profile, while the same changes 
would result in linear changes in the HSE mass estimate, 
as can be seen by examining Equation 1101 

Adopting the Na varro, Frenk, and White profile 
([Navarro et al.lll997l hereafter NEW) to describe the to- 
tal mass distribution (i.e., baryonic 4- dark matter dis- 
tribution), the total density is radially distributed as 



Ptot(r) 



(r/P,)(l-f r/P,)2' 
and the total mass within r is 



(11) 



= AnpoRl [ln(l + r/P,) - (1 + Rs/r)-'] . 

(12) 

Here po a-nd are respectively the normalization and 
scale radius of the NEW profile. The use of an NEW 
profile is empirically motivated by simulations of dark 
matter halos, and we note that other mass profiles could 
be assumed and may in fact provide better alternatives. 
More recent theoretical studies, for instance, have indi- 
cated that the presence of baryons can sign ificantly mod- 
ify the mass distributio n of dark matter ([Gnedin et al.l 
[200l[mldd"erall[200l . 

The gas mass Mgas and total mass A/tot can be related 
by defining the gas fraction /gas = Mgas(r)/Mtot(r), 
which could be a function of Mtot, r, z, and cluster 
merger histo ry. While recently shown to be poor app rox- 
imation (e.g. IVikhlinin et al.|[200l iPratt et al.ll2010l) . we 
make the simplifying assumption that fgas{r) is a con- 
stant. Detailed measurements of /gas('") typically require 
high significance X-ray data, which are often insufficient 
or entirely lacking for clusters discovered via the SZE. 
The assumption of constant /gas(f) implies that 



Pgas('") = /gasPtot(r-)- 



(13) 



We assess the impact of this assumption in Section [4] 
Using Equations (TTJ [121 & [HI to solve for the gravita- 



3 



tional potential energy, we find 
Ug{r) =(47rpoi??)'G/gas 



[2{1 + Rjry 



'■ ln{l + r'/R,) 
(l + r'/i?,)2 



dr' 



(14) 



where we have used the fact that the differential ele- 
ment of gravitational energy for a spherical shell of gas 
with density pgas{r) is dUg{r) = —GMtotdm/r, where 
the mass of the gas shell is dm ~ Air p^i^s{f')r'^ dr . 

Combining EquationslTlfc lT^ through the virial relation 
(Equation in]) , we have 



(1 + V/^e) 
167r2G/ga. 

Rs 



^rUeC 



'■ln(l 



2(l + iZ,/r)2 



r'/Rs 



(l + r'/i?,)' 



■dr' 



(15) 

Here the —Airr^ Pc{r) term in the brackets on the left 
hand side of the equation is due to the surface pressure 
correction in Equation |9j 

Through the above relation one can find the best fit 
NEW profile parameters po and Rs for any observation- 
ally constrained Ysp\i{r). We apply this method to inter- 
ferometric SZE data in Section |3| 

3. APPLICATION TO OBSERVATIONS WITH THE 
SUNYAEV-ZEL'DOVICH ARRAY 

We test here the application of Equation |T5] to the 
SZE observations of the three clusters presented in 
IMroczkowski et all (120091 hereafter M09), and compare 
our results with the independent mass determinations 
presented in M09. These three clusters span a wide 
range in redshift and dynamical state. A1835, at z = 
0.25, is a relaxed, cool-core cluster (e.g., iPeterson et al.l 
1200 It) . A1914, at z = 0.17, shows evidence of be- 
ing d isturbed, with a hot subcluster near the cluster 
core (|Maughan et al.l l2008l ). CL J1226.9-H3332 (z = 
0.89) appears somewhat rel axed given its high redshift 
(|Maughan et al.ll2004 120071 ). but recent high- resolution 
SZE o bservations with MU STANG have indicated oth- 
erwise (|Korngut et al.ll2010D . 

M09 derived mass estimates for these clusters using 
three data fitting methods. The first method relied on 
SZE observations and X-ray surface brightness data, but 
ignored the X-ray spectroscopic data. Instead, a density 
model was fit to the X-ray surface brightness data si- 
multaneously with a pressure profile fit to the SZE data. 
Temperature information used in fitting the X-ray sur- 
face brightness data were derived from these density and 
pressure profiles, assuming the ideal gas law. The den- 
sity model used in this method was a core-c ut simpli- 
fication of that used in iVikhlinin et al.l (|2006l hereafter 
V06), to which we refer as the "Simplified V06 Model" 
(SVM). The pressure profile used in this method is an 
analytic parameter ization of th e clust er radial pressure 
profile proposed bv INagai et al.l (|2007l hereafter N07), 



Peir) 
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This profile has the form of a generalized NEW profile, 
and was fit with the slopes fixed at the best fit values 



found in N07, which are (a,6,c) = (0.9,5.0,0.4)0 We 
refer to this X-ray-|-SZE method, which does not use or 
require X-ray spectroscopic data, as the "N07+SVM" 
method. 

The second mass estimation method presented in 
M09 was an independent, X-ray only analysis per- 
formed by Ben Maugha n following t he methods out- 
lined in iMaughan etall (|2007l 120081) . This state of 
the art method relies on deep Chandra X-ray observa- 
tions, fitting both the spectroscopic and surface bright- 
ness data with the full density and temperature param- 
eterizations in V06. The X-ray only analyses of each 
cluster in M09 were published in thr ee separate pa- 
ers: CL J1226. 9-^3332 was publ ished in IMaughan et ed] 
2007D . A1914 was published in IMaughan et all (12008^ 
and Al835 was published in M09. 

The third method fit the SZE and X-ray data jointly, 
but relied on the assumption of isothermality and was 
included in M09 only for comparison with earlier works. 
We ignore this method here, noting that its results were 
consistent at small radii but increasingly discrepant at 
large radii. 

With the density profile and the temperature or pres- 
sure profile from the above methods in hand, the equa- 
tion of hydrostatic equilibrium (Equation llOp was used 
to solve for each cluster's total mass. The results of the 
N07+SVM and X-ray only analyses are reproduced in 
Table m 

We compare the results of the M09 mass analyses with 
those from the SZE-only method presented here. We as- 
sume the same ACDM cosmology used in M09 {flM = 
0.3, Qa = 0.7, and h = 0.7). As in M09, we adopt the 
N07 profile ( Equation HEl w i th (a, &,c) = (0.9,5.0,0.4)) 
and that from lArnaud et all (|2010l . hereafter AlO, Equa- 
tiondiwith (a,6,c) = (1.0510,5.4905,0.3081)). 

The SZE data used here were taken with the Sunyaev- 
Zel'dovich Array (SZA). Briefly, the SZA is an eight el- 
ement compact array built to image the SZE in clus- 
ters through observa tions at 30 and 90 GHz (see, e.g., 
iMuchovej et al.ll2007| ). At 30 GHz, the SZA is sensitive 
to radial scales l'~6' over a 10'.6 diameter field of view. 
At 90 GHz, the SZA measures 20"-120" radial scales 
over a 3'. 5 diameter field of view. These observations are 
naturally fit in u, w-space (Fourier space) using a Markov 
Chain Monte Carlo (MCMC) process, as discussed in 
M09. The trial model is computed in the plane of the 
sky by integrating Equation [3l using Equation [16] to de- 
scribe the pressure. We present here the results assuming 
both the N07 and AlO parameterizations. The sky-plane 
model is then Fourier transformed for direct comparison 
with the interferometric data. This has the advantage 
that the likelihood of the model fit is computed in a ba- 
sis where error bars arc Gaussian. The full details of this 
method can be found in M09. 

The resulting fit N07 and AlO pressure profiles are used 
to solve for the radial profile of Ysph{r) (Equation [5]). For 
each accepted link in the MCMC, we fit the function de- 

^ The parameters published in N07 are (a, b, c) = (1.3, 4.3, 0.7), 
but the combination (a,b,c) = (0.9,5.0,0.4) was later found to 
provide a better fit. The planned erratum to N07 is yet to be 
published, and the corrected parameterization first appeared in 
M09. However, the corrected parameterization we use here has 
come to be known as the "Nagai 07 Profile Pressure," and we 
adhere to this terminology. 
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TABLE 1 

Kjph AND Aftot FOR Each Model Tested, Computed within r25oo and rsoo Assuming Constant 

/gas = 0.13 



Cluster Name 
Model Fit 



''2500 
(Mpc) 



'5^sph(''2500) 
(lO-^Mpc^) 



Mtot(r2500) 

(IQI^Mq) 



''500 

(Mpc) 



i'sph(''500) 

(lO-^Mpc^) 



Mtot{r50o) 

(IQi^Mo) 



A1835 

N07 (this work) 
AlO (this work) 
N07+SVM (M09) 
Maughan (M09) 

CL J1226+3332.9 
N07 (this work) 
AlO (this work) 
N07+SVM (M09) 
Maughan at al. (2007 ) 

A1914 

N07 (this work) 
AlO (this work) 
N07+SVM (M09) 
Maughan et al. (2008 ) 



0.63 
0.63 
0.68 
0.66 



+0.01 
0.01 
+0.01 
-0.01 
+0.02 
-0.02 
+0.02 



0.39 



+0.01 

-0.01 



40+°-°l 
'-'■^'J-0.01 

0.4ll-«l 

0.60 
0.59 



+0.02 
0.02 
+0.02 
0.02 

^•"'-0.03 
,+0.02 



0.63" 



0.02 



7, 
7, 

8.25 
7.88 

3.34 
3.54 



+0.52 
-0.50 
+0.48 
-0.50 
+0.81 
0.78 
+0.49 
0.72 



+0.28 

-0.28 
+0.28 



-0.29 
o t:f5+0.36 
■^■"^"-0.36 

5 04+°-^^ 



4.59 
4.37 



+0.75 
0.65 
+0.69 
0.57 



6 29+1"^ 
"■^^-0.82 

-1+0.37 



5.69" 



0.38 



4.58'' 



j+0.18 
-0.19 

^■°*-0.19 
5 RA+0-5S 



5.30 



0.54 
+0.53 
0.72 



r, or:+0.15 
^■'^''-0.16 

r, c-n+0.14 

^■'-''-'-0.15 

Z.D/_Q 27 
o 4. +0.30 
■^■^^-0.26 



3.59 
3.49 



+0.30 
0.30 
+0.31 
0.29 



4 07+0.89 
-0.72 
,+0.43 



4.31" 



0.33 



1.45 
1.44 
1.44 
1.42 

0.94 
0.94 



+0.04 
0.04 
+0.04 
0.04 
+0.11 
0.10 
0.07 
0.05 



+0.02 
0.02 

+0.02 



-0.02 

0.98t°;J? 
0.89tro2 



1.27 



+0.07 
0.07 



1 90+0.07 
^■'^■'-0.06 

1.25 



1.29 



0.10 

+0.07 
0.06 



17.66 
16.21 



+2.87 
2.36 
+2.47 
2.07 



+3.00 
1(.00_2 70 



17.41 



+ 1.61 



9.43 



+0.89 

-0.85 



Q 17+0.88 
^■^'-0.83 
,+1.58 



1.29 



9.71] 



10 59+0-69 



7.77 



+2.29 
1.70 



6.67+1;^^ 
10.78±i;°| 



11.24 
10.87 



+ 1.03 
-0.94 
+0.99 



ii.ool^-g 

10.68l?'^f 



6.49 
6.42 



+0.34 
-0.34 
+0.36 
-0.36 



0,7+2.50 

'■■"-1.57 
5 4Q+0.46 



6.88" 



-1.26 
-1.07 



6.62tM° 

3+1.29 
-1.00 



7.49"* 



scribed by Equation [TSl assuming a constant /gas =0.13 
and fie = 1-17. The resulting mass profiles, computed 
using Equation [121 are used to find ta and Mtot (''A ) , 
which are respectively the radius within which the aver- 
age density is A times greater than the critical density 
of the universe at that redshift, and the total mass con- 
tained within that radius. As in MOO, wc report ta, 
Mtot(rA), and ysph(rA) for A = [2500,500], with statis- 
tical error bars, in Table [TJ Wc discuss systematics along 
with our conclusions in Sectional 

4. CONCLUSIONS 

We conclude that this method is remarkably consistent 

— given the simplifying assumptions required to derive 
total mass using the virial relation and SZE data alone 

— with the X-ray only and X-ray-f-SZE mass determina- 
tion methods in M09 (sec Table [ij . The assumption of 
constant /gas = 0.13 has perhaps the largest systematic 
impact on the derived values of -Mtot (''a) and ta for over- 
density A. The radial mass profile -MtotCT) is oc /gas~^^^, 
as can be seen by examining the relation between the 
NEW parameter po and /gas in Equation [TH However, 
any change in Mtot (f ) affects ta and therefore Mtot (''A ) , 
so the systematic change in the mass at fixed overdensity 
is larger than a simple rescaling by the inverse square root 
of the ratio of the correct to the assumed /gas- Eitting 
the same data with an assumed /gas — 0.11, for example, 
increases Mtot (''a) by an average of 12% (rather than 
the 9% change in the profile Mtot(''))- 

The assumption that fie = 1-17, by contrast, can be ex- 
pected to have a much smaller impact on the mass deter- 
mination method presented here. Eor typical abundance 
gradients due to metal enrichment, fic varies on the '^1% 
level, which changes the (1 -I- l//io) factor in Equation 
[TSl on the ^ 0.5% level. Large systematic deviations in 
metallicity therefore affect the fit Mtot('') at the ^ 0.25% 
level. The assumption of a single, constant metallicity is 



also common in X-ray studies of high redshift clusters, 
where the limited number of X-ray counts is insufficient 
to constrain more than a single spectroscopic bin. On 
the other hand, helium sedimentation in the absence of 
magnetic fields, in a cluster undisturbed fo r 3 Gyrs, could 
increase He in the core re gion by ~ 5% (iPeng fc Nagail 
I2009| ). Using the results of lPeng fc Nagail |2009(). we note 
that the sedimentation of helium has little effect on the 
average fic or on fJ-dr) at large radiiH The expectation 
is that mergers and r nagnetic fields will b oth suppress 
helium sedimentation ([Peng fc Nagail [2009( ). 

Another potential source of bias is due to uncertain- 
ties in the calib r ation of the SZE data. As discussed in 
IMuchovei et al.l ([20071) . the absolute calibration of SZA 
data is known to better than 10%, and the variation from 
observation to observation in amplitude of a flux cali- 
brator (in this case Mars) is < 5%. Calibration errors 
would result in scalar systematic errors in the fit pres- 
sure profile and have a linear impact on Vsph (Equation 
[5]). Examining the relation between Vsph and the NEW 
parameter po in Equation I15[ we can see errors in the 
derived Mtot(?') oc ysph('')^^^, and will impact the mass 
at fixed overdensity -Mtot(?'A) on the < 5% level. 

As noted in M09, SZA 30 GHz observations are sen- 
sitive to radial angular scales ^1-6', so the largest scale 
measured is ~ r^oo for A1835, and is ~ 0.87-500 for A1914. 
The values we report should therefore be treated as ex- 

* In fact, the effect of helium sedimentation is greater for X-ray 
surface brightness data than for SZE data. The redistribution of 
helium nuclei into the core greatly increases the X-ray emissivity, 
impacting both the charge and number of ions, and the number 
of electrons, the product of which determines the bremsstrahlung- 
dominated X-ray emission. The intensity of the SZE is only im- 
pacted by a factor proportional to the increase in the number of 
electrons. Therefore, the mass determinations from X-ray data can 
be expected to be more biased by the effects of heliu m sedimenta- 
tion than those based on the SZE. For good reviews, see lMarkevitchI 
ll2007l) : IFeng fc Nagail | |2009|') . 
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trapolations of the fit, and depend on the assumed N07 
or AlO pressure profile. 

The surface pressure term in the virial relation (Equa- 
tion [9l) was not included in the original work, and can 
be of the same order as the thermal energy of the clus- 
ter gas. While this term systematically reduces the mass 
inferred through this method, the amount the estimate 
changed from the value neglecting surface pressure was 
typically within the error bars. Overall, the inclusion of 
this term was found to improve the agreement between 
the SZE-only inferred mass and the estimates from X-ray 
and SZE-hX-ray. 

Other systematics the SZE-only mass estimation 
method could suffer are common to X-ray mass deter- 
minations that rely on a cluster's fit radial temperature 
and density profiles and assume thermal HSE to derive 
mass. First, HSE is most readily applied by assuming 
spherical symmetry. Second, while the gas may be viri- 
alized within the potential and supported predominantly 
by thermal pressure, there is an expectation that '^10%- 
20% of th e total pressure is due to turbulent mo tions in 
the ICM (|Lau et all [200l iBattaglia etUIIMol ). This 
is equivalent to including kinetic energy, in addition to 
thermal, in the virial relation. The interesting question 
arises as to whether adopting the virial relation, instead 
of just HSE, could be used to make X-ray mass determi- 
nations more robust. 



The broader implications of this work are that by solv- 
ing for the radial thermal energy profile one can esti- 
mate cluster mass from any SZE observation that can 
constrain that cluster's radial pressure profile. Radial 
profiles have already been fit to clusters observed by 
ACT, SPT, and n ow Planck (iMenanteau et a 1 120101 : 
iPlagge et all I2OIOI: iPlanck Collaboration et all l2011aD . 
and this method may prove particularly useful for provid- 
ing initial mass estimates for the clusters they discover. 
A future work, using SZA observations of a more com- 
plete sample of clusters, will compare total mass derived 
in this way with X-ray and lensing mass estimates. We 
will also consider how to extend this method for the case 
where the gas fraction varies with cluster radius. 
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